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Abstract. Elegant and mathematically rigorous methods of the quantum in- 
verse theory are difficult to put into practice because there is always some 
lack of needful input information. In this situation, one may try to construct 
a reference potential, whose spectral characteristics would be in a reasonable 
agreement with the available data of the system's properties. Since the refer- 
ence potential is fixed, it is always possible to calculate all its spectral charac- 
teristics, including phase shift for scattering states and Jost function, the main 
key to solve the inverse problem. Thereafter, one can calculate a Bargmann 
potential whose Jost function differs from the initial one only by a rational 
factor. This way it is possible, at least in principle, to construct a more re- 
liable potential for the system. The model system investigated in this paper 
is diatomic xenon molecule in ground electronic state. Its reference potential 
is built up of several smoothly joined Morse type components, which enables 
to solve the related energy eigenvalue problem exactly. Moreover, the phase 
shift can also be calculated in part analytically, and the Jost function can be 
acertained very accurately in the whole range of positive energies. Full energy 
dependence of the phase shift has been determined and its excellent agreement 
with the Levinson theorem demonstrated. In addition, asymptotically exact 
analytic formulas for the phase shift and the Jost function, independent of 
each other, are obtained and their physical background elucidated. 



1. Introduction 

Strict mathematical criteria for the unique solution of the inverse problem for 
the one-dimensional Schrodinger equation have been formulated more than 50 years 
ago. The history of this prestigious research area goes back to the dawn of quantum 
mechanics 1 , and the decisive breakthrough has been achieved thanks to the im- 
portant contributions by Borg 003, Levinson Bargmann [HJ|7], Gel'fand and 
Levitan [Hj, Jost and Kohn jHHEj, Marchenko [III El, Krein OGl, and others 
(see, e.g., |15| for a thorough overview). As a result of fruitful brainstorming which 
culminated in early fifties it has been established that the quantum mechanical 
inverse problem can be solved if one manages to fix the so-called spectral func- 
tion jl lj . In the case of confining one-dimensional quantum systems it means that 
the interaction potential can be uniquely determined if and only if the following 
complete set of information is available: 

• Full energy spectrum of the bound states E n (n = 0, 1, N). 

• Full energy dependence (from to oo) of the phase shift 6(E) for the scat- 
tering states E > 0. 
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• n additional real parameters C n (n — 0, 1, N) related to relevant bound 
states that uniquely fix their normalization. 

In principle, all these data can be obtained experimentally, but unfortunately, 
this is almost unachievable in practice. The real situation is even more hopeless, 
because in addition to the deficit of input information one inevitably faces very 
serious computational-technical problems. Thus, in spite of the whole mathematical 
beauty of the theory, one comes to a regrettable conclusion that rigorous solution 
of the quantum-mechanical inverse problem is a tremendously difficult task. For 
this reason most methods of deducing potentials from the available experimental 
data are based on some simplifying preconditions, therefore being inaccurate from 
the rigorous quantum mechanical point of view. Nevertheless, such methods may 
prove quite useful. For example, semiclassical approaches introduced many years 
ago by Rydberg [HUH], Klein E3> Rees and Dunham [2U E21 are still very 
popular in spectroscopy of diatomic molecules, and these concepts are constantly 
improved and developed |23l I24j . 

Is it possible to apply rigorous methods of the inverse quantum theory for prac- 
tical purposes? In view of the principle difficulties mentioned above, one has to be 
cautious in answering this question. In this paper we restrict ourselves to simple 
one-dimensional quantum systems, and the analysis proceeds from an idea that for 
any system of this kind one can build up a reasonable reference potential based on 
the available experimental data. As the reference potential is known, it is always 
possible to calculate all its discrete energy eigenvalues E n , their norming constants 
C n , and the phase shifts 6(E) for the scattering states. Therefore, in this artificial 
way one gets the full set of input information needed to uniquely solve the quan- 
tum mechanical inverse problem. Of course, the described approach is tautological: 
there is no need to regain a potential which is already known by definition. Never- 
theless, such an approach is not meaningless, as it gives good zeroth approximations 
to the important spectral characteristics, such as Jost function and spectral den- 
sity (the terms to be specified below) . One can interpret the reference potential as 
only an initial guess to the real potential. Although the calculated quantities E n , 
C n and 6(E) do not exactly match the actual values for the real system, they are 
still expected to be quite close to them. For a given reference potential one can 
calculate its Bargmann potential whose Jost function differs from the initial one 
only by a rational factor ^S]. By a suitable choice of this factor one can take a 
more adequate account of the experimental data. For example, one can replace the 
calculated discrete energy eigenvalues E n related to the reference potential with 
their actually observed values. Consequently, at least in some sense the new po- 
tential would be more realistic than the initial reference potential. There is also 
another motivation for the described "inverse" approach to the inverse problem. 
Namely, through direct practical experience one can essentially increase his knowl- 
edge of how to overcome serious computational-technical difficulties when applying 
rigorous methods of the inverse quantum theory. 

The one-dimensional inverse theory can be applied to diatomic molecules, since 
a two-particle problem can be always reduced to a one-particle problem in a spher- 
ically symmetric field. On this basis, we are going to examine the inverse problem 
for diatomic xenon molecule in its ground electronic state. Several reports of the 
research are planned. Methods of solution of the integral equations that enable to 
uniquely ascertain the potential are discussed and illustrated in the next paper of 
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this series [lfi| . while in this paper the emphasis is put on explaining the details 
of the basic concepts related to the reference potential approach. In Section 2 we 
briefly describe how the reference potential for the model system has been con- 
structed, and how its discrete energy eigenvalues have been calculated. In Section 
3 we describe the details of calculating the phase shift for the scattering states, and 
demonstrate full agreement with Levinson theorem g|. Section 4 aims to explain 
the important role of the Jost function in the quantum inverse theory. In particular, 
a detailed analysis of the asymptotic behavior of the Jost function is given. Finally, 
a brief conclusion is given in Section 5. 

2. Exactly solvable reference potential for Xe 2 

In Fig. 1 one can see the reference potential constructed for the Xe 2 molecule. 
The same curve is depicted in both graphs, but very different energy scales are 
used. Throughout this paper only the rotationless case is analyzed, i.e., the ro- 
tational quantum number J = 0. According to the starting idea of the approach 
the only criterion for the choice of the reference potential is its agreement with 
the available experimental data. Therefore, we will not pay too much attention 
to various mathematical nuances and simply assume that the reference potential 
should be smooth and integrable in the whole physical domain. In addition, we 
try to construct a reference potential whose analytic form is as simple as possi- 
ble. A good choice for this purpose, as explained in detail elsewhere [251 I26| . is a 
multi-component potential composed of smoothly joined Morse-type pieces 

(1) V(r) =V k + D k [exp(-a fc (r - r k )) -lf,r€ (0, oo), 

where Vk, D kl a k and r k are some constants (not definitely positive), and the 
subscript k corresponds to different components smoothly joined at some suitably 
chosen boundary points -Xfc.fi. The reference potential shown in Fig. 1 consists 
of only three components [k = 0, 1, 2), the most internal of them (k — 0) being a 
so-called pseudo-Morse potential. It means that the tiny potential well correspond- 
ing to this component (if taken separately) is just of the limit depth to entirely 
lose the discrete energy spectrum. Consequently, D = ^ 2 0!q/(8to) (to being the 
reduced mass of the pair of atoms), so that only three independent parameters 
remain for this component. The central component (k = 1) is an ordinary Morse 
potential, while the most external one (k = 2) is a "reversed" Morse potential with 
the parameter _D 2 being negative. By introducing a "reversed" component one 
artificially creates a small potential hump in the long-distance range. This might 
seem unphysical and unjustified, but the point is that the height of this artificial 
hump approaches zero as the parameter r 2 approaches infinity. Therefore, taking 
a sufficiently large r 2 , the hump becomes almost insignificant, while the analytic 
treatment remains simple and flexible. All parameters of the reference potential 
can be easily determined, if one requires continuity of the potential and its first 
derivative at the boundary points X\ and A 2 (also shown in Fig. 1). These pa- 
rameters as well as the calculated discrete energy eigenvalues E n (n = 0, 1, 23) 
are given in Table 1. 

The essence of the described construction is that the energy eigenvalue problem 
for the reference potential can be easily solved by solely analytic means to any 
desired accuracy [22] ■ Moreover, as we demonstrate in the next Section, the major 
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part of the energy dependence of the phase shift can also be ascertained analytically, 
which is a great advantage compared with applying numerical methods. 

Most spectroscopic applications are related to the distance region shown in the 
lower graph of Fig. 1. For the inverse quantum theory, however, it is important 
to accurately reproduce the potential near the zero point r = 0, which is almost 
meaningless for spectroscopic applications. In this context one cannot ignore the 
fact that according to Eq. (1) the reference potential is finite at r = (see the 
upper graph of Fig. 1). Actually it means that the potential "jumps" to infinity at 
zero point, and this might also seem unphysical and unjustified. However, one has 
to bear in mind that the behavior of the real potential near r = is unknown and 
remains unknown. It does not matter so much how we describe the potential in this 
region, in so far as spectroscopic applications are of our main interest. As we see in 
the next section, description in terms of a pseudo-Morse potential is mathematically 
simple and elegant, and this, too, can be taken as a motivation for the approach to 
be used. 



3. Phase shift and Levinson theorem 

Next we are going to calculate the phase shift 8(E) for the full range of scattering 
states E € (0, oo). To this end, one can use a long-known method first introduced by 
Morse and Allis |57] (let us remind that only the rotationless case J = is examined 
here). It is based on solution of the following equation: 

-/2mV(r) . 2 



(2) S'(r, k) = - - w sin 2 [kr + S(r, k)] , 6(0, k) = 0, 

nk 



where k = — . The phase shift is then determined as S(k) = lim 5(r, k). The 

described method is universal but rather time-consuming, because for any energy 
from to oo one has to perform an integration from to oo. Fortunately, at this 
point we can take advantage of the special analytic form of the reference potential. 
As mentioned, the region r S (0, Xi) is approximated by a pseudo-Morse potential. 
It means that one immediately gets two linearly independent analytic solutions of 
the corresponding Schrodinger equation |251 I26| 

(3) ^\r) = exp(-y /2)2/r ^0,2*^0 + l',yo) 

* 2) (r) =e X p(y /2)(-yo) iM *(iMo,2iMo + l;-yo), re(0,X 1 ), 

where ^ = 1/2 VC^ ~ Vo)/D - 1 > 0, y = exp [-a (r - r )], and *(i// , 2i/io + 
1; t/o) is a particular solution of the confluent hypergeometric equation introduced 
by Tricomi (see |2H| for details). If /ig << yo, this function can be evaluated from 
the asymptotic series 



(4) #(i/i ,2iMo + l;yo) = %^ ^ 



(iya)n(-ino) 
n\(-y ) n 



where (a) n = T(a + n)/T(a) = a(a + l)(a + 2). ..(a + n — 1) is the Pochhammer 

symbol, and N must not be too large. Thus 

(5) 

Mo , Mo(Mo + l 2 ) Mo(Mo + 1 2 )(Mo 



^(r) = exp(-y /2) 



l!yo 2!y 2 i\yl 
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and analogously 

(6) - exp(W2) [l + 4- + + iM+^Oilg + , 

The phase shift is related to regular solutions of the Schrodinger equation, which 

means that the physically correct linear combination of "J/j, (r) and (r) should 
vanish as r — > 0, i.e., 

y (r)=N 1 ^\r)+N 2 ^ ) (r), 



where 



M§ + Mo(Mo + I 2 ) 



-i , Mo , ms(mo + 1 



l!yo(0) 2!y 2 (0) 

(2) 

From Eq. (7) one can see that the particular solution <3> (r) practically does 
not contribute to the regular solution at distances sufficiently far from zero point, 
but still deep inside the classically forbidden region. Indeed, in this case 

1 | Mo^ + l 2 ) 

, n , N 2 ^\r) . . . , n ., l!y o (0) 2!y 2 (0) 

( 8 ) m =-exp[y (r)-t/o(0)] x 



+ l!y o (0) + 2! % 2 (0) 

Mo , Mo(Mo + l 2 ) , 



x *fl> = ~ exp M r) - ^o(O)] x 

Mo | MovMo "t" 1 



/'o 



i!yo(0) 



l!y (r) 2!j/g( r ) 
[exp(a r) - 1] + ... 



i!yo(0) 



- exp {- exp(a r ) [1 - exp(-a r)]} 



is an extremely small quantity. 

There is still a lot more profit to gain from the pseudo-Morse approximation to 
calculate the phase shift. As we just proved, the physically correct solution (apart 
from normalization) in a wide energy range (practically up to E ss V(0)) reduces 
to the particular solution (r). Using some well-known relations from the theory 
of confluent hypergeometric functions J211 IZ@ this solution can be rewritten 

(9) *^(r) = exp(-y /2)j, w, vi/(^ 0; 2j Mo + l;ifo) = 



exp(-y /2) 



r(-2z/i ) iM0 r(2jMo) - V(W . „. L1 v 

— — — -y^ $(«/i ,2i/x + l yo) + ^ — r-% $ -*Mo, -2«Mo + 1 J/o) 



a a; a(a + 1) x 

where the symbols $(a, c; i = H H ) — - + ... denote the well-known 

c 1! c(c + 1) 2! 

confluent hypergeometric functions. 

As is seen, Eq. (9) represents a sum of two complex conjugates. Therefore, 

*o ( r ) ~ ^o(yo) cos [B (yo) -<po- ao^ r] , 
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where the functions Aq and Bq can be determined from a series |25| : 

A (yo)e-°W = 1 - + ^'*>* U - -^-) 

i^ + l/2 (z/x + 1/2)1! \ «/ i o + 3/2y 

(10) + . (yo/4) 4 A yo/4 



(i/io + V 2 ) («Mo + 3/2) 2! V »Mo + 5/2 

A good point is that the phase parameter (which, of course, is not the actual phase 

T(2i M o)" 



shift) tpo = a ^or - arg 



can be calculated exactly. Indeed, using a 



r(z/i ) 

Legendre formula for doubling a gamma function's argument |28| 

r(2z) = ^pr(z)r(z + i/2), 

and another useful formula [2H| 
(11) 

OO 

argr(i/xo+l/2) = Q ln(l + 4p§) - 1 - ln2j ~j (cothi - ^ e~* sin(2 Mo i)^, 

o 

one gets the following result: 
(12) 

OC 

¥?o=Mo ^a r + l-ln2-iln(l + 4^)^ +i ^ (cotht - ^ e _t sia(2/i <) j) 

o 

where the integral can be conveniently evaluated |25| 

oo T 

(13) / = y (cotM - e~* sin(2 Mo i) j = J sm(w^)f(t)dt, T = JL, 



= cotht-i _ cothft + T)-^ cothft + 2T)-^ _ 

f 6 i + T 6 t + 2T 

Another equivalent formula for this quantity reads 26 



/ x (-l) n - 1 2 2n B„ , fc /2n-lV 



,2fc+l 
0) 



Here B n denotes the n-th order Bernoulli number. 

Correct linear combinations of solutions of the Schrodinger equations for analyt- 
ically different pieces of the reference potential are uniquely fixed by the continuity 
requirements of the wave function and its first derivative at the boundary points 
X\ and X2- Since the parameter tpo can be calculated exactly, the phase shift 
can also be ascertained with the help of solely analytic means in the energy range 
< E < V(0). At high energies (E > V(0)), however, the analytic approach grad- 
ually fails, since the particular solution ' (r) in Eq. (3) cannot be ignored any 
more. To determine the phase shift in this region one has to solve Eq. (2). For- 
tunately, there is no need to perform integration over the whole physical domain 
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r € (0,oo), because the solution in the long-distance range r > X 2 already has the 
"right" analytic form [2H] 

(15) * 2 (r) = 2A 2 (y 2 ) cos [B 2 {y 2 ) + <p 2 - kr] = 2A 2 (y 2 ) sin [S(r, k) + kr] , 

from which the phase shift can be easily obtained. Here y 2 (r) = 2a 2 exp [— a 2 (r — r 2 )] 
a 2 = y / 2mD 2 / (ha 2 ) , S(r,k) is the solution of Eq. (2), and the complex function 
A 2 (y 2 ) exp [iB 2 (y 2 )} = exp (-iy 2 /2) $ [i(k/a 2 - a 2 ) + 1/2, 2ik/a 2 + 1; iy 2 \ . Since 
A 2 (y 2 ) — > and B 2 (y 2 ) — > 1 as r — ► 0, the phase shift 

(16) 5(k) = 8(X 2 ,k) + B 2 [y 2 (X 2 )] , 

i.e., one only needs to integrate Eq. (2) until the boundary point X 2 . 
Now, let us recall an important relation known as Levinson theorem 0] 

(17) (5(0) - 6(00) = mr, 

which correlates the energy dependence of the phase shift with the number of bound 
states. As can be seen in Fig. 2, a really good agreement with the Levinson the- 
orem can be obtained, but only if the phase shift is calculated up to very high 
energies (note that the energy scale in Fig. 2 is logarithmic, and it involves 20 
orders of magnitude!). Moreover, the phase shift has to be calculated with suf- 
ficiently high precision throughout the whole energy range, otherwise there is no 
chance to accurately ascertain other important spectral characteristics, such as 
Jost function. Unfortunately, the higher the energy goes, the more complicated 
and time-consuming the numerical integration of Eq. (2) becomes. How could we 
bridge over this troublesome technical difficulty? In such situation, one may recall 
some general principles, and this is indeed helpful to complete calculations of the 
phase shift. 

Let us express the phase shift as formal result of integration of Eq. (2): 
(18) 



1 



JT V(r)dr 



8(k) = lim 5(r, k) = — / V(r) {1 - cos [2kr + 26 (r, k)]} dr 



2Ck J 1 K ' 2Ck 

OO 

1 



.,, .cos (2kr) {V(r) cos [25 (r, fc)]} dr - — — / sin (2kr) Mr) sin [25(r, k)}} dr. 
2Ck J 2Gk J 



Here a special denotation has been introduced (and will be used henceforward) for 

the constant C = - — that often appears in formulas. Now, let us call for help 
2m 

from the famous Riemann-Lebegue theorem (see, e.g., [30]): if a function F(r) is 
integrable in an interval r € (a, b), then 

b b 

(19) J cos (Ar) F(r)dr ->• 0, J sin (Ar) F(r)dr -> as A -> 00. 

a a 

From this one immediately concludes that the last two integrals in Eq. (18) will 
vanish as k — > 00. However, we can integrate by parts and apply the Riemann- 
Lebegue theorem to the resulting integrals. This procedure can be repeated as 
many times as needed, and as a result, one comes to a rather general asymptotic 
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formula for the phase shift (it can easily proved that only odd powers of k appear 
in this series) 

(20) S {k) = % + % + % + ..., fc^oc, 
where 

/ V(r)dr ( oo 

(21) ai = - ° , a 3 = -gij | J [V(r)} 2 dr + CV'(0) 



«5 = - 32^3 I / i V (r)} 3 + dj [V'(r)] 2 dr + 4CV(fl)V'(fl) - C 2 V"'(0) 
I o o 

Thus one can easily calculate the coefficients a%, a 3 , 05, etc., and this is just what 
is needed to ascertain the whole energy dependence of the phase shift. Fig. 3 
demonstrates how well Eq. (20) fits with the results of direct numerical integration 
of the phase equation. 

One cannot so easily find any direct illustrations to the Levinson theorem from 
the literature (at least the author of this paper has not found them), and Fig. 2, 
which is an illustration of this kind, could therefore be of more general interest 
than merely an attachment to a particular model potential. The curve shown in 
this figure has been calculated with at least 8 significant digits, and it is in full 
agreement with all relevant general physical considerations. This concerns not only 
the Levinson theorem and the asymptotic behavior at large energies, but also the 
low energy part of the energy dependence (see the left-side inset of Fig. 2), exactly 
corresponding to the well-known formula 

(22) S{k) = fin - arctan(fci? ), k -> 

(i?o being scattering length), which can be found in most handbooks on quantum 
mechanics. 



4. JOST FUNCTION AND INVERSE PROBLEM 

Having calculated the phase shift, we have come to a situation from which an 
ideal inverse problem study would start. In other words, we are now provided with 
the full set of information needed to uniquely solve the inverse problem. In our 
case this would mean that we simply regain the reference potential from which 
we started. This, of course, is not our main goal. As explained in Section 1, we 
are interested in getting realistic zeroth approximations to the important spectral 
functions which then could be used to improve the initial model potential. The 
latter step, however, is planned as a subject for a forthcoming publication. In this 
paper we only try to make sure that the described scheme is reliable, and to this 
end the next step is to accurately calculate the Jost function, the main spectral 
characteristic, which contains the most part of information needed to solve the 
inverse problem. A thorough overview of all useful properties of the Jost function 
is given in Ref. |15|. For the treatment here, the most important point is that the 
Jost function creates a simple link between physical and regular solutions of the 
scattering states. The physical solution (J = 0) reads (cf. with Eq. (15)) 

(23) *(r, k) w expiS(k)sm[kr + 6(h)] , r — > 00, 
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while the regular solution ip(r, k) is defined by a condition 
(24) 



lim ^ = 1. 



i — >o r 

It can be shown that these two solutions of the Schrodinger equation are proportional 

k 



(25) 

where 
(26) 



\E»(r, k) 



F(k) 



ip(r, k), 



F(k) = |F(fc)|exp [-iS(k)] 



is the Jost function we are talking about. For further treatment we have to calculate 
the modulus of this function ^5] 



A 



(27) 



\F(E)\ = l[(l-E n /E)exp 



n=0 



1 f 5{E')dE 

ir J E' -E 
o 



E e (0,oo) 



Here E n are the discrete energy eigenvalues and the symbol P points at the principal 
value of the integral. As the phase shift and the bound states are known, the 
calculations are relatively trivial, but they have to be carried out very accurately. 
Next we can fix the spectral density 



(28) 



dp{E) _ J n^VElFjE^ 2 ,E>0, 
dE ~ 1 J2C n S(E-E n ), E<0. 



Here C n are the norming constants for the relevant bound states. Note that these 
quantities are related to regular solutions, wherefore their ascertainment is not so 
easy task as one might think. 

Now we have come close to the real solution schemes of the inverse problem. For 
example, the Gelfand-Levitan method [5j is based on the integral equation 



(29) 

whose kernel reads 
(30) 



K(r, r') + G(r, r') + J K(r, s)G(s, r')ds = 0, 
b 



, , sin (fcr) ■ sin (fcr') 
G(r,r)= I - da, 



fc 2 



and the quantity da = dp(E) dE. Here - is free particle's spectral 

density (related to the potential V(r) = 0), and therefore, 



2R3/2 

dp(E) — d(— ), E >0 
da = { J7T 

dp(E), E < 0. 



(31) 

Eq. (30) can be rewritten [T5] 

oc 

2 r (j 

(32) G(r, r') = — / sin (fcr) ■ sin (fcr') g(k)dk + > — ^ sinh (7„r) sinh ( r y n r l ) , 

* J „ 4 7 2 



10 



MATTI SELG 



where 7^ = — ^ m ^ n an( j ^ f unc tion 

in h 2 



(33) 



9(k) = 



\F(k)\< 



- 1. 



If one is able to solve Eq. (29), the potential can be determined from the relation 



(34) 



V(r) 



Eqs. (29) to (34) explicitly demonstrate great importance of the Jost function in 
the inverse quantum theory. The energy dependence of the Jost function's modulus 
is shown in Fig. 4. One can see that at small energies this quantity achieves 
extremely large values. Near the "critical" energy E = V{Q) the curve rapidly 
turns from nearly vertical to nearly horizontal, and at still higher energies it slowly 
approaches the limit value 1^(^)1 ■— > 1, "breaking free" from the potential field 
(note that for a free particle, F(k) = 1 independent of energy). 

The kernel of the Gelfand-Levitan equation is essentially determined by the func- 
tion g(k) given by Eq. (33), which means that this function has to be ascertained 
very accurately, in order to solve Eq. (29) and calculate the potential according 
to Eq. (34). As can be seen in Fig. 5, g(k) noticeably differs from unity only at 
k > fc = y/V(0)/C, and g(k) — > as k — > 00. Naturally, as shown in the inset, 
there is no break of derivative in the "critical" region. Since the asymptotic be- 
havior of the function g(k) essentially determines the shape of the potential near 
the zero point r = (see the end of this section), it makes sense to analyze this 
behavior in more details. We have already ascertained the asymptotic expression 
for the phase shift (see Eq. (20)), and this can be used to immediately get the 
asymptotic formulas for both the Jost function's modulus and the function g(k). 
For example, taking 



(35) 



1 m/.M a 2 a 4 a e 
ln|F(fc)| = - + F + - + ..., 



k 



00 



(this time only even powers of k appear in the series, as can be easily proved), the 
coefficients read 

(36) 



(Z4 



:(aik a 



2 aifc| 

7T l 3 



03 

k a 



«5 

3fc3 



•) + 



f 



azk a 



05 

k a 



...) + 



C 2 



J S(E')dE' -J2 E n 



J 5(E')E'dE'- l -Y,^y 



2 ai kl 
a 6 = --(—=— 

7T 5 



+ + a 5ka 



•••) + c3 



6(E')(EfdE'--J2(E n f 



where E a — Ck\ is an arbitrary energy value in the range where the asymptotic 
approximation Eq. (20) can be used. 

Eq. (36) may look nice but it is a bit inconvenient in practice. Fortunately, a 
straightforward approach exists, which enables to ascertain the coefficients in Eq. 
(35) more accurately and much more easily without any direct reference to the phase 
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shift. The approach in question is based on the following integral representation 
for the Jost function 



oo 

(37) F(k) = 1 + i J e* kr V(rMk : 



r)dr, 



where the regular solution ip(k, r) can be calculated with the help of a well-known 

sin kf 

iteration method. Namely, taking tp^(k,r) — — ; — , and 

k 

r 

M (h ^ = I f sinfc(r-rQ ■ („_!) , , 



(38) V w (k,r) = -I y - >-V{ r ')^-V{k,r')dr', n = l,2,.., 

o 

the desired solution reads 

oo 

(39) i P (k,r) = Y / 9 {n) (k,r). 

n=0 

In Eq. (38) one can use integration by parts to calculate step-by-step the terms 
(/^^(fc, r), (p( 2 \k,r), <pw(k,r), and their contributions to the Jost function. Let 
us see, how to ascertain the correct asymptotic formula for the function g(k) = 

x — 1 as k — > oo. which would include the terms until ~ l//c 4 . Within this 

\F(k)\ 2 . 

approximation 
(40) 

r 

V il} ( k >r) = -^W(r) + S ^[V(0) + V(r)] + [ V > (r) - V (0)] , W ( r ) = J V{r')dr> , 



y (2) (fc>r) = - Sinfc ^y )]2 - [V(0)W(r) + V(r)W(r) + U(r)] , U(r)= \[V{r')fdr', 

(3) _ coskr-[W(r)f 
V {k ' r >- 48CW 

(all higher order terms can be ignored). Thereafter, using integration by parts in 
Eq. (37), one comes to the following formulas for the real and imaginary parts of 
the Jost function: 

(41) 

V(0) W 2 V"(0) 1 r o 1 

RGF ^ = 1 + S - 8C*k 2 - iW + 32C^ I 5 ™? - 2V ^ W ) - 

' [V(0)W 2 + 2UW] U 



32C 3 fc 4 L w J 384C* 4 fc 4 ' 



W V'(0) 1 w 3 
(42) Im F(k) = — + -A-i + -_ [V(0)W + U\- 



2Ck 8Ck 3 8C 2 k 3 J A8C 3 k 3 ' 

OO OO 

where W = J V{r')dr' and U = J [V(r')f dr' . Quite surprisingly, when calculating 

o o 
\F(k)\ 2 = [ReF(fc)] 2 + [ImF(k)] 2 , all troublesome terms will cancel out, resulting 
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in nice expressions for both quantities of interest: 

V(0) 3[V(0)} 2 V"(0) 



(43) \F(k)\< 



2Ck 2 8C* 2 fc 4 8C*fc 4 ' 



(44) m .m + w^m. 

Comparing Eqs. (35) and (44), one finds that the coefficients 

F(0) 2 [V(0)f - CV"(0) 

(45) a 2 = ^,a 4 = _ 

do not depend on the phase shift, but are directly related to the potential and its 
second derivative at the zero point, both these quantities being finite according to 
the starting idea of the approach. Thus we have described all properties of the 
Jost function in the frame of the proposed approach, and are now prepared to start 
calculation of the potential itself. 



5. Conclusion 

In this paper we proposed as if an " inversed" approach to the quantum mechan- 
ical inverse problem. Starting from a known reference potential, one can calculate 
the important spectral characteristics of the system, the phase shift and the Jost 
function, which are almost unattainable in a real experiment. On one hand, the 
reference potential has to be realistic enough to be used as a zeroth approximation 
to the real potential. Only in this case there is a chance to construct a Bargmann 
potential that would be even more realistic, for example, whose discrete energy 
levels would exactly fit with the actually observed ones. On the other hand, we 
suggest to choose a reference potential that would be exactly solvable, in the sense 
that its energy eigenvalue problem can be solved to any desired accuracy with the 
help of solely analytic means. To this end, as demonstrated in Sections 2 and 3, 
a multi-component potential composed of smoothly joined Morse-type pieces is es- 
pecially suitable. In particular, we would like to highlight the usefulness of the 
pseudo-Morse approximation for the small distances (and high energies) region of 
the potential. In this paper we only used a single pseudo-Morse component, which 
stretches until the zero point r = 0. As demonstrated elsewhere [2E1GE], adding 
more pseudo-Morse components does not bring along any serious problems, so one 
can include just as many components of this type as he considers reasonable. 

The special analytic form of the reference potential enables to ascertain the phase 
shift analytically up to the energies E j§ V(Q) (note that an arbitrarily large value 
for V(0) can be taken, if one introduces a sufficient number of pseudo-Morse com- 
ponents). As the asymptotic behavior of the phase shift can be determined from 
general physical considerations, it is possible to directly demonstrate full agreement 
with the famous Levinson theorem 0], as well as to ascertain the full energy de- 
pendence of the Jost function. Provided with this important input information, 
one can attack the real computational-technical problem, trying to solve the main 
integral equation in the frame of the Gelfand-Levitan [H] , Marchenko or Krein 
[T4*j method. This will be the subject of the next paper in this series |16) . 
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reference potential approach to the inverse problem: i 15 
Figure captions 

Fig. 1. Three-component model potential for the system (Xe2 in ground electronic 
state) investigated. Note that the same potential curve is depicted in both 
graphs (the lower one starts where the upper one ends), but essentially dif- 
ferent energy scales are used for them. All components have the well-known 
analytic form of the Morse potential, but the ordinary Morse approxima- 
tion is used only in the central range r G [Ai,X 2 ] (sec the explanations 
in Section 2). The parameters of the components as well as the calculated 
discrete energy levels (24 in total) are given in Table 1. 

Fig. 2. Direct demonstration of the Levinson theorem (8(0) — S(oo) = nw) for the 
model system studied. As needed, 6(0) = 24tt, since the system has 24 
bound states. At E — 3.146294 meV, the phase shift passes a zero, and 
then remains negative, very slowly approaching the limit (6(oo) = 0) as 
E — > oo. The left side- inset shows the nearly linear energy dependence as 
E — > 0, in full agreement with Eq. (22). In the right-side inset one can see 
that the phase curve has an inflection point near E = V(0). 

Fig. 3. A demonstration of the exellent agreement between the numerically calcu- 
lated phase shifts and the general asymptotic formula, Eq. (20), as E —* oo. 

Fig. 4. Demonstration of the results of calculating the Jost function (in fact, its 
modulus) for the reference potential. Again, as in the case of Fig. 1, 
different parts of the same curve are shown in the graphs and their insets, 
but rather different scales are used (note that in some cases the scale is 
logarithmic). The information gathered into the upper graph would be 
useful, if one's aim is to construct another potential whose Jost function 
only slightly differs from the initial one. The lower graph demonstrates 
how abruptly the curve turns from practically vertical to nearly horizontal 
in the vicinity of the "critical" energy value E = V(0), while this dramatic 
change is absolutely invisible in the scale used for the ordinate axis in the 
upper graph. The asymptotic formula Eq. (35) is used to calculate the 
Jost function for k > k a = 75000 A (which corresponds to about 1.8*10 5 
eV) with coefficients given by Eq. (45). 

Fig. 5. Visualization of the characteristic function g(k) (given gy Eq. (33)), which 
determines the kernel of the Gelfand-Levitan integral equation. As can 
be seen, g(k) — 1 (with very high precision) if k J$ k — y/V(0)/C (see 
Fig. 4 to undestand the reason for this), then starts to decrease, gradually 
"breaking free" from the potential field, and g(k) — > as k — > oo. In spite 
of the seeming simplicity of the energy dependence, the whole curve has 
to be (and has been) calculated very accurately up to very high energies, 
in order to accurately ascertain the potential. The high-energy part of the 
curve (k > k a ) has been calculated according to Eq. (44). 
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Table 1 . Parameters of the three-component reference potential 
for Xe2 (first column) and its discrete energy eigenvalues E n . 
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